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We study analytically and numerically the winding of directed polymers 

of lenght t around each other or around a rod. Unconfined polymers in pure 
media have exponentially decaying winding angle distributions, the decay con- 
stant depending on whether the interaction is repulsive or neutral, but not 
on microscopic details. In the presence of a chiral asymmetry, the exponen- 
tial tails become non universal. In all these cases the mean winding angle is 
proportional to Int. When the polymer is confined to a finite region around 
the winding center, e.g. due to an attractive interaction, the winding angle 
distribution is Gaussian, with a variance proportional to t. We also examine 
the windings of polymers in random systems. Our results suggest that ran- 
domness reduces entanglements, leading to a narrow (Gaussian) distribution 
with a mean winding angle of the order of v / m~7. 

I. INTRODUCTION AND SUMMARY 

The topological constraints produced by the windings of polymers Q strongly affect the 
dynamics of polymer solutions. As a consequence of polymer entanglement, the viscosity of 
a solution of polymers above the overlap concentration is many orders of magnitude higher 
than the viscosity of the solvent. An analytical treatment of these topological constraints 
is extremely difficult, and theoretical efforts therefore focus on the limit of high polymer 
concentrations, where effective medium theories and the tube model successfully describe 
several aspects of the dynamics of the polymer solution 0, or on the limit of only one or 
two polymers, where the different possible configurations can be studied explicitely ||. 

In this article, we take the latter approach, focussing on the winding of a directed polymer 



(DP) around a rod, or of two DPs around each other, as shown in Fig. [I]. DPs have a preferred 
direction f, and their configuration can be described by the function {f(r)}, with r G [0, t], 
where r — (£1,22) is the coordinate in the plane perpendicular to the preferred direction. 
Going to relative coordinates r^ir) — r*i(r), the winding of two DPs around each other can 
be mapped to the winding of a single DP around a rod (see section [TT]). Since DPs cannot 
have knots, their main topological constraints are windings. 

Although less common than flexible polymers, DPs are a good model of several semi- 
flexible and rigid polymers. Example are biological macromolecules such as DNA or liquid 
crystals composed of stacks of disk-shaped molecules. These polymers are aligned parallel 
to each other when their concentration is sufficiently large, forming crystalline and liquid 
crystalline phases (for a review on statistical mechanics of DPs see e.g. [f|]). Isolated DPs 
can be realized by embedding a long polymer in a nematic solvent ||. Another important 
class of directed "polymers" are magnetic flux lines in high-T c superconductors that are ori- 
ented parallel to the direction of the external magnetic field. Due to the high temperature 
in the system and the weak coupling between different layers in the superconductor, thermal 
fluctuations of the flux lines are considerable, leading to entanglements H0] and windings 
around columnar pins ||. 

In order to calculate the winding angle distribution of DPs, we map them onto two- 
dimensional walks {^(t)}, where the arc length r plays the role of the time coordinate. 
The winding angle distribution depends on the interaction between the polymer and the 
winding center, and the properties of the embedding medium. In the following three sec- 
tions, we discuss three different classes of winding angle distributions, for each of which the 
scaling variable is a different combination of the winding angle 9 and the polymer length 
t. In section [TI], we consider DPs in an infinitely large pure medium. These polymers can 
be mapped on ideal random walks, the mean horizontal distance (|f(r) — ^*(0)|) from the 
starting point increasing with the square root of r. The number of returns to the winding 
center is proportional to lnt for such a random walk. We will see below that because of 
the finite return probability to the winding center even for large times, the winding angle 
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distribution depends on properties of the winding center. We will find different winding 
angle distributions depending on whether the interaction between the winding center and 
the polymer is a hard-core repulsion or is neutral, and whether the winding center shows 
a chiral asymmetry. However, the winding angle distribution does not depend on micro- 
scopic details like the shape of the winding center or a possible underlying lattice structure. 
The probability distribution for the winding angle depends in the limit t — > oo only on 
the combination x = 20/ ln(t) of the winding angle and the length of the walk, and not 
on each variable separately. For large \x\, all three mentioned winding angle distributions 
decay exponentially in \x\. The scaling variable proportional to 9/ ln(i) can be explained as 
follows: After time t, the walker has a typical distance r(t) oc \/i from the starting point, 
which is chosen to be close to the winding center. Assuming that r(t) is the only relevant 
length scale, dimensional arguments, combined with the Markovian property, suggest that 
dr/d9 = rf(9). The rotational invariance of the system implies that f{9) must be a con- 
stant, i.e. the increase in winding angle cannot depend on the number of windings or angular 
position. Hence, 

d9 oc y oc j = d(lnt) , (1) 

leading to a scaling variable proportional to 9/ Int. 

If the interaction between the polymer and the winding center is attractive, the polymer 
can be bound to the winding center, and its transverse wandering is limited. Polymers can 
also be confined by a finite container or by neighboring polymers. In all these cases, polymer 
segements of length At that are small compared to the total length t, but large compared to 
the length needed to make a winding, have identical winding angle distributions. Applying 
the central limit theorem, we conclude that the total winding angle distribution is a Gaussian 
with a scaling variable 9 2 /t. This situation will be discussed in section |T|. 

Finally, we discuss in section [TV] certain DPs that cannot be described by ideal random 
walks, since they do not satisfy a Markov property. Polymers that are embedded in a random 
medium, e.g. a gel, have an energy that depends on the polymer configuration. Similarly, 
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the energy of magnetic flux lines in high-T c superconductors with point defects depends 
on their configuration. The mean distance from the starting point for these polymers in 
random media increases faster than for ideal walks, since the line searches for low-energy 
configurations. The number of returns to the winding center remains finite in the limit of 
infinite length. Consequently, the interaction of the line with the winding center does not 
affect the winding angle distribution, as long as it is not strong enough to bind the polymer. 
We will see that the winding angle distribution of polymers in random media is a Gaussian 
with a variance proportional to hit, leading to a scaling variable x oc 9/ J\a{t). This means 
that the pinning to randomness decreases the mean winding angle from the order of In t to 
the order of Vint. Interestingly, this winding angle distribution is similar to the one for two- 
dimensional self-avoiding random walks |J. The following scaling argument explains why 
the winding angle distribution in both situations is Gaussian with a variance proportional 
to hit: Starting from the origin divide the walk into segments of 1, 2, • • • , 2 n rs t/2 steps. 
Since the a th segment is at a distance of roughly 2 au from the center (with v — 3/4) and has 
a characteristic size of the same order, it is reasonable to assume that each segment spans a 
random angle 9 a of order one. Under the mild assumption that the sum 9 = J2a=i @a satisfies 
the central limit theorem, we then conclude that 9 is Gaussian distributed with a variance 
proportional to n oc In t. Since this argument relies on the irrelevance of the winding center, 
it cannot be applied to the distributions in section 0. 

Many results of this article have been reported previously in Ref. ||10|| . They point out 
the rich behavior already present in the simplest of problems involving topological defects. 
Properties of the winding center, interactions, various types of randomness are all potentially 
relevant, leading to different universal distribution functions. The concluding section [V| of 
this article gives an outlook on possible further universality classes and on the winding of 
non-direceted polymers. 
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II. WINDING ANGLE DISTRIBUTIONS IN AN INFINITE HOMOGENEOUS 



MEDIUM 

In this section, we study winding angle distributions of DPs in infinite homogeneous 
media, all characterized by a scaling variable x = 26 / hit and exponential tails. We keep 
the initial point of the polymers fixed, but otherwise allow them to move freely. The precise 
form of the winding angle distribution depends on the interaction with the winding center. 
In subsection [11 A|, we consider two DPs with hard-core repulsion, or, equivalently, one DP 



winding around a repulsive rod, leading to the distribution in Eq. (f[4|) . For neutral winding 
centers, the corresponding winding angle distribution given in Eq. ( |17|) has a decay constant 
that is smaller by one half (subsection |II B| ). These two distributions occur under fairly 



general conditions (see subsection [II C|) . However, when the symmetry with respect to the 
sign of the winding angle is broken, new (asymmetric) distributions occur (subsection |II D| ) , 
with the decay constants of the exponential tails depending on the degree of chirality. 



A. Winding in the presence of hard-core repulsion 



1. Mapping to a random walk with absorbing boundary conditions 



The energy of a given configuration of two DPs ri(r) and r 2 (r) of length t is given by 



E[f 1 (r),f 2 (r)]= fair 
Jo 



dr 



dr 



(2) 



The potential V(r) has a hard-core, V(r) = oo for r < a, and V(r) = for r > a. The first 
two terms are the elastic energies of the polymers, where the parameter c is related to their 
stiffness. Introducing the relative coordinate f ' = r\ — r 2 and the center-of-mass coordinate 
R= [ri + f 2 )/2, Eq. (g) becomes 



E 



R(r),r(r) 



dr 



c ( dr\ I dR . 
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The partition function for the two polymers is 



Z = Jv [R(t)] V [f(r)] exp [-E [r(t), r(r)] /k B T} 



(3) 



where the integral is taken over all possible configurations [-R(t)] and [r(r)]. The expression 
V[r{r)] denotes a path integral and is the continuum limit of nr=i(/ dr(Ti)), k B is the 
Boltzmann constant, and T is the temperature. 

As long as we are only interested in quantities related to the relative coordinate, like the 
winding angle, we can integrate out the center-of-mass variations, and focus on the partition 
function for the relative coordinate alone, i.e., 



T> [r(r)] exp 



dr 



/k B T . 



(4) 



This is identical to the partition function for a single DP winding around a rod. Due to 
the hard-core repulsion, all configurations where the polymer and the rod penetrate each 
other, do not contribute to the partition function (V = oo), while V = for all other 
configurations. 

A two-dimensional random walk can be described by the Langevin-equation 



(5) 



where rj is a stochastic force with zero mean ((rj(t)) = 0) and the correlation function 
(rj(t)ff(t')) = 2D5(t — t'). The probability distribution of ff is a Gaussian, i.e. 

P [r){t)\ oc exp [-D (ff(t)f . 

With Eq. (H), we find that the probability for a given trajectory [f(r)] of the random walk 
is proportional to 

exp 



When all walks that enter a region of radius a around the origin get absorbed, the probability 
that the random walk has a trajectory [r(r)] is identical to the probability that the above 
DP has the cofiguration [r*(r)] (compare to Eq. (^), with D = c/ksT). This correspon- 
dence between DPs with hard-core repulsion and random walks with absorbing boundary 



conditions was first pointed out by Rudnick and Hu [11 
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2. Conformal mapping of the random walk 



Since we are interested in the winding angle of the random walk, it is convenient to 
perform a transformation such that the winding angle becomes one of the coordinates. To 
this purpose, we represent the walk r{t) = (xi(t),x 2 (t)) by the complex number 

z(t) = xi(t) + ix 2 (t) . 

The time evolution of each random walker satisfies 

dz = fj(t)dt, (6) 

where rj (t) is now complex, with 

( V (t) V *(t')) = 2D5(t-t>). (7) 
We now introduce the new variable £ by the transformation 

C(t) = lnz(t) =\nr(t)+i9(t), (8) 



where p = lnr = In y x\ + x\. Since d( = r](t)dt/ z(t), the stochastic motion of the walker in 
the new complex plane is highly correlated to its location, i.e., the walk is no longer random. 
This feature can be removed by defining a new time variable 

dt 

dt = WW (9) 

for each walker, which leads to 

d( = n(t)dt, with /i(t) = z*(t)r}(t). (10) 

Since 

(//(£)//(?)) = 2D\z(t)\ 2 5(t - t') = 2D5{t-t'), (11) 

the evolution of ((t) is that of a random walk. Under the transformation in Eq. (|3|), the 
absorbing disc in the plane z maps onto an absorbing wall in the plane ( (see Fig. H). 



For simplicity we choose the initial condition ((t = t = 0) = 0, i.e. the original walker 
starts out at z — 1. We also set the diffusion constant to D — 1/2, so that the mean 
square distance over which the walk moves during a time t is (r 2 (t)) = t. Consequently, the 
probability that r(t) is within an interval [y/nt^ 1 ' 6 ^ 2 , y/Tit^ 6 ^ 2 } around its mean value of 
y/nt, is 



and approaches unity in the limit t — > oo. The effect of the absorbing disc on this probability 
can be neglected in the limit t — > oo, since the disc becomes smaller when viewed from 
larger distances. In this limit, the distance r from the starting point z = 1 is identical to the 
distance from the origin, and p(t, e) is identical to the probability that ((£) is in the interval 
[0.5(1 — e) In t, 0.5(1 + e)lnt]. So the endpoints of all walks (except for an infenitesimal 
fraction) that take a time t in the original plane map within a strip of width elnt in the 
(■-plane, as indicated in Fig. [| If we shrink the complex plane ( by a factor of (lnt)/2, 
the walker is within a distance e of the line with real value of unity. Thus, all walks of 
length t in the z-plane are mapped on walks that end at the line with real value of unity, 
without having gone beyond (see Fig. 0). Since there is a separate transformation t(t) for 
each walker, walks of the same length t map on walks of different length t. (To be precise, 
we also have to shrink the time scale t when shrinking the (-plane, but for simplicity we 
denote the new time again by t.) 



The imaginary (or vertical) coordinate x in the rescaled (-plane is related to the winding 
angle by x = 26(t)/\nt. In order to obtain the winding angle distribution, we have to 
determine the vertical position of a random walk starting at the origin, at the moment 
when it reaches for the first time the wall at distance one from the vertical axis, without 




(12) 



3. Calculation of the winding angle distribution 



S 



going beyond the absorbing wall at distance 2| lna|/lnt on the opposite side. Since we are 
interested in a walk only up to the moment when it reaches the right-hand wall, we can 
consider this wall also as absorbing. Since walks of length t in the original plane map on 
walks of different length t in the new plane, we need the probability that the walk is absorbed 
at this wall before time t. 

We formulate this problem more generally and determine the probability P a ^(y,t) that 
a one-dimensional random walk starting at y between two absorbing points a and (3 at time 
is absorbed at the point /3 before time t. Since for sufficiently small At, the walker is only 
a short distance Ay from its starting point, we have 



r°° 1 
P a Av,t) = / d{Ay)-== 
Jo V27rAt 



cxp 



2At 



x 



P aJ} (y + Ay, i - At) + P a>p (y - Ay, t - At) 



Expanding the above equation to the order of At indicates that P a ^(y, t) satisfies a diffusion 
equation. The appropriate boundary conditions are P a ^(a,t) = and P a ^{j3,t) = 1 with 
the initial value P a> p(y, 0) = 0, resulting in fP2" 



P a ,p(y,t) 



y-a 



(3 — a 
x exp 



r> oo 



-l)^ 1 (ixviy-a) 
sm 1 



(5 — a 



1 



7T1S 



t 



The probability that the walk is absorbed at the right-hand boundary during the time 



interval 



t,t + dt 



, is dtdfP aj p(y,t). 

Note, however, that / °° dtdiP a ,p(y,t) = (y — ot)/(f3 — a), i.e. equal to the total fraction 
of particles absorbed at the right-hand boundary (inversely proportional to the separations 
from the boundaries). To calculate the winding angle distribution Pa{x), we need the fraction 
of these walks absorbed between t and t + di, equal to — a)/(y — a))d^P. Hence (with 
a — 2 In a/ In t, /3 = 1 and y — 0) 



p A (x) 



00 d ~ 1 - a dP aA (0, t) exp (-x 2 /It 
o -a dt ^2~Kt 
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poo x 



-1 



v 2vrt «(1 - a 

1 / 7TV 



TXV ( Trva 

sin 



a 



E 

i/=i 



x exp 
1) 



a 



2 ~ x 2 

2 V 1 - a ) ~ 2t 

in/a \ 

sin [ exp 

1 — ay 



irv\x\ 



1-a 



The last step is achieved by first performing a Fourier transform with respect to x, followed 
by integrating over t, and finally inverting the Fourier transform. (Alternatively, the t 
integration can be performed by the saddle point method.) In the limit of large t, the 
variable a is very small, and we can replace the sine-function by its argument. 
Taking the sum over u, we find 

exp [irx/(l — a)] 



p A {x) 



[1 - a) {exp [ttx/(1 - a)] + 1} 



2 ' 



(13) 



Changing the variable from x to 



x 



26 



x 



(1-a) ln(t/a 2 )' 
and noting that pa(x)cLx = pA(x)dx, leads from Eq. ([Dp to 

( 26 \ 



7T 



Pa x 



(14) 



ln(t/a 2 ); 4cosh 2 (vrx/2) ' 

The above distribution, which is exact in the limit t — > oo, has an exponential decay 
for large x, as first derived in Ref. [IT]. The complete form of Eq. (O) was first given 
in Ref. f[3"|] , however, without derivation. The analogy to random walkers in the plane (, 
confined by the two walls, provides simple physical justifications for the behavior of the 
winding angle. In the presence of both walls, the diffusing particle is confined to a strip, 
and loses any memory of its starting position at long times. The probability that a particle 
that has already traveled a distance 6 in the vertical direction proceeds a further distance 
dO without hitting either wall is thus independent of 6, leading to the exponential decay. 
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4- Comparison to the winding angle distribution around a point center 



The method described in this section was used earlier to derive the winding angle distribu- 
tion for Brownian motion around a point center [13]. The resulting probability distribution 
for the winding angle in this case is [13] 

leading to an infinite mean winding angle. Since there is no confining wall on the left-hand 
side, the particles may diffuse arbitrarily far in that direction, making it less probable to 
hit the wall on the right hand side. In the original z-plane, the walker takes no time at all 
to make an infinitely small winding around the point center. This is clearly an unphysical 
feature, since real winding centers are finite and since real random walks (or polymers) need 
a finite time (a finite length segment) to make a winding. We therefore do not consider this 
situation any further. 

Exercise: Derive Equation ftTh]), repeating the calculation of this section, but with no 
absorbing wall (see Ref. JTZj/j. 



B. Winding of directed polymers around neutral winding centers 

Instead of having first a rod and then inserting the DP into the system, we can also 
first have a free configuration of a DP, and then insert a rod into it. If the polymer cannot 
relax to its thermal equilibrium distribution after insertion of the rod, e.g. because its 
end are fixed or because its configuration is frozen, the resulting winding angle distribution 
will be different from that in the previous subsection. No configuration of the polymer 
is forbidden, but those configurations that interfere with the rod become deformed. The 
degree of deformation may depend on the diameter of the rod, but the winding angle does 
not. Alternatively, we could consider a winding center that has no interaction at all with the 
polymer, e.g. a light beam, or some structural defect in the solvent that is not felt by the 
polymer. In this case we would find the same winding angle distribution as in the case of a 
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rod that deforms the polymer. In the language of a random walk, this situation corresponds 
to having a disc that reflects all walks that hit it. The walks that would go through the disc 
thus become deformed, but are not removed from the statistical ensemble. 

We can obtain the winding angle distribution by repeating the calculations of the previous 
subsection, but replacing the absorbing boundary condition P a ^(a, t) — with the reflecting 
condition dP a ^{y,t) /dy \ y=a = 0, leading to 

p / r\ , 2 " 1 . f n(u + l/2)((3-y) \ 
n^ Q v+l/2 \ (3 -a J 



x exp 



1 /tt(z/ + 1/2)\ 2 



(16) 



2\ (3 -a J 

There is thus no current leaving the system at point a, and walkers which hit the winding 
center are reflected. We then find the winding angle distribution 

VR ^ = 2coshU/2) ' (17) 

where again x = 26/ ln(i/a 2 ), and the limit t — > oo has been taken. For large x, where 
the walk has lost the memory of its initial distance from both walls, this probability decays 
exponentially as exp [—ttx/2], i.e. exactly half as fast as for absorbing boundary conditions. 
A random walk confined between an absorbing and a reflecting wall that have a distance 
one can be mapped to a random walk confined between two absorbing walls at distance two. 
After rescaling the wall distance and the x-coordinate by two, this explains the factor 1/2 
between the decay constants in the tails of the distributions in Eqs. ( |14D and (fl7|). 



Exercise: Derive Equation (11) (see Ref. fdW) 



C. Universality of the winding angle distribution 



The winding angle distribution in Eq. (|ITD, which we derived in the previous subsection 
for Brownian motion around a reflecting disc, was obtained previously by several authors in 



different contexts: Belisle []TB| calculated the winding angle distribution for a random walk 



on a two-dimensional lattice around a point that is different from any lattice site, and for 
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a random walk with steps of finite size taken in arbitrary directions, around a point in the 
two-dimensional plane, obtaining in both cases Eq. fllTD- The same result was obtained by 
Pitman and Yor [Tl for the distribution of "big windings" of Brownian motion around two 



point-like winding centers. Comtet, Desbois, and Monthus (18| divided the two-dimensional 
plane into three concentric sections and determined the contribution of each section to the 
windig angle for Brownian motion around a point, finding Eq. (|17|) for the contribution of 
the outer section. This universality seems surprising, since one might expect that the main 
increase in winding angle occurs when the walk is close to the winding center, where details 
like the lattice symmetry, and shape and size of the winding center determine how much time 
it takes to make one winding. However, a careful look at Fig. |2| reveals that this is not the 
case: The main increase in winding angle does not occur when the walker is within a small 
distance from the left-hand wall. Since all distances have been scaled by l/ln(t), a small 
distance from the left-hand wall corresponds to a large distance (of the order of In t) from 
the winding center. Therefore, almost all windings are made far from the winding center, 
where microscopic details do not matter. The properties that do affect the winding angle 
distribution are conservation laws (absorbing or reflecting boundary conditions), symmetries 
(with respect to the sign of the angle - see the following subsection), singularities (as for 
the winding of Brownian motion around a point center), and interactions (self-avoidance) 
or randomness (see section |TV|) . 

To further test this unversality hypothesis, we determined numerically the winding angle 
distribution for a random walk on a lattice with reflecting and absorbing boundary condi- 
tions. Reflecting boundary conditions are realized by choosing a winding center different 
from the vertices of the lattice, and thus never crossed by the walker (this is exactly the 
situation treated analytically in Ref. ||16|| ). On the other hand, to model absorbing boundary 
conditions, the winding center is chosen as one of the lattice sites (say the origin), but no 
walk is allowed to go through this point. 

The winding angle distributions are most readily obtained using a transfer matrix method 
which calculates the number of all walks with given winding angle and given endpoint after 
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t steps, from the same information after t — 1 steps. The winding center is at (0.5, 0.5) for 
reflecting boundary conditions, and at the origin for absorbing boundary conditions. The 
walker starts at (1,0), and the winding angle is increased or decreased by 2n every time 
it crosses the positive branch of the xi-axis. Due to limitations in computer memory, we 
applied a cutoff in system size and winding angle for times t > 120, making sure that the 
results are not affected by this approximation. The largest times used, t = 9728, required 
approximately 3 days to run on a Silicon Graphics Indy Workstation. 

Figures ||] and |] show the results for the two cases. The asymptotic exponential tails 
predicted by theory can clearly be seen; deviations from the theoretical curve for smaller 
values of the scaling variable x = 28/ ln(2t) are due to the slow convergence to the asymptotic 
limit. Since the scaling variable depends logarithmically on time, the asymptotic limit is 
reached only for large Int. Note that the only free parameter in fitting to the analytical 
form is the characteristic time scale appearing inside the logarithm. With t measured in 
units of single steps on the lattice, we found that a factor 2 in the scaling variable provides 
the best fit. In the limit t — > oo, different scales of t give of course the same asymptotic 
winding angle distribution. 

Exercise: Perform the numerical calculations mentioned in this section. Study also the 
case of absorbing and reflecting winding centers that comprise several lattice points. How 
does the size of the winding center affect the convergence towards the asymptotic winding 
angle distribution? 

We also studied the winding of a DP proceeding along the diagonal of a cubic lattice 
in three dimensions (see Fig. ||). The polymer starts at (1,0,0), and at each step increases 
one of its three coordinates by 1. We determined the winding angle distribution around the 
diagonal (1, 1, 1)- direction, excluding from the walk all points that are on this diagonal (a 
repulsive columnar defect, corresponding to the case of an absorbing winding center). The 
excluded points lie on the origin when the polymer is projected in a plane perpendicular to 
the diagonal. In this plane, the polymer proceeds along the bonds of a triangular lattice, 
alternating between the three different sublattices. A cutoff of 243 in system size was imposed 
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for the transfer matrix calculations. The winding angle distribution p(x = 28/ hit) is shown 
in Fig. ||for different times. As for the square lattice, an exponential tail with decay constant 
of 7r can be seen. Our numerical results, as well as the analytical considerations, thus indicate 
clearly that the winding angle distributions for reflecting and absorbing boundary conditions 
are universal and do not depend on microscopic details. 

Due to the special properties of directed paths along the diagonal of the cube, the case 
of reflecting boundary conditions leads to an asymmetry between windings in positive and 
negative directions. This is because it takes only three steps to make the smallest possible 
winding in one direction, but six steps in the opposite direction. This situation is discussed 
in detail in the following subsection. 

D. Winding centers with chiral asymmetry 

So far, we only considered situations that are symmetric with respect to the angles ±9. 
For directed paths on certain lattices, however, this symmetry is broken. A directed walk 
that proceeds at each step along +21, +22, or +X3 direction on a cubic lattice can be mapped 
on a random walker on a two-dimensional triangular lattice as indicated in figure 0(a). Each 
bond can be crossed only in one direction, and the winding center for reflecting boundary 
conditions must be different from the vertices of the lattice. It is apparent from this figure 
that the random walker can go around the center in 3 steps in one angular direction, but in 
no less than 6 steps in the other direction. 

An alternative description is obtained by examining the position of the walker after 
every three time steps. The resulting coarse-grained random walk takes place on a regular 
triangular lattice, but now the walker has a finite probability of 3 x 2/3 3 = 2/9 of staying at 
the same site. If this site is one of the three points next to the winding center, the winding 
angle is increased by 2tt in one of the six possible configurations that return to the site 
after three steps. In other words, the walker has a finite probability of having its winding 
angle increased in the proximity of the center. The amount of this biased increase in angle 
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depends on the structure of the lattice and will be different for other directed lattices. An 
equivalent physical situation occurs for Brownian motion around a rotating winding center, 
e.g. a rotating reflecting disc which does not set the surrounding gas or liquid into motion 
(see Fig. §). 

Since this angular symmetry breaking is already present in the above simple example 
of a directed walk, it is quite likely to occur in more realistic physical systems, such as 
with screw dislocations in the underlying medium. The winding angle distribution for two 
chiral polymers f| should also show an angular asymmetry. We thus use the term chiral 
winding center to indicate that each time the polymer comes close to the center, it finds it 
easier to wind around in one direction as opposed to the other. (Of course, to respect the 
reflecting boundary conditions, there must be no additional interaction with the winding 
center. The case of additional attractive or repulsive interaction with the winding center 
will be discussed later.) 

After mapping to the rescaled £-plane introduced in Sec. [IT], the above situations can be 
modeled by a downward moving, reflecting wall on the vertical axis. Each time a random 
walker hits this wall, its vertical position x = 26/ \nt is changed by a small amount 2 AO/ Int. 
Let us now determine the net shift Ax in x due to the motion of the wall for a walker that 
survives for a time t in the rescaled £-plane, before it is absorbed at the right-hand wall 
(recall that t is the time in the original system, while t refers to the time in the rescaled 

plane, after the conformal mapping). 

To obtain the full solution, it is necessary to solve the two-dimensional diffusion equation 
with moving boundary conditions. Since we are mainly interested in the exponential tails 
of the winding angle distribution, we restrict our analysis to the limit of large times t, and 
determine the shift in x due to encounters with the reflecting moving wall in this limit. A 
Brownian walker that has survived for a sufficiently long time t forgets its initial horizontal 
position. The mean number of encounters with the reflecting wall, and consequently the 
shift Ax in x due to the motion of the wall, is then expected to be simply proportional to 
the considered time interval. Applying the central limit theorem in the limit t — > oo, the 

16 



probability distribution of Ax is given by 

Pa(Ax) = 



27T(3 2 t 



exp 



(Ax - aif 
2(5 2 t 



(18) 



The parameters a and (3 are related to the velocity of the wall (chirality of the defect) by 
a oc (3 oc v . Presumably Eq. can be obtained directly from properties of random walks, 
providing the exact coefficients for these proportionalities. 
The tail of the winding angle distribution is then given by 

(x — Ax) 2 



p c R (x)oc! dt I d(Ax) dPa i ( l ),t) i-^ xp 



dt 



oc / dtexp -(7r 2 /4)t/2 



'2wt 
1 



27rt(l + (3 2 ) 



exp 



2t 

(x - at) 2 
~2t{l + [3 2 ) 



exp 



ax 



\x\ 



7T Z 



1 + f3 2 1 + (3 2 



l a 2 + — (l + (3 2 ) 



(19) 



valid for large \x\. The second factor on the right-hand side of the first line of Eq. ([19]) is 
the probability that the random walker is absorbed at the right-hand wall at time t (see 
Eq. (|i6|), and is for large t dominated by the slowest mode iy = 0). The second factor 
is the probability that the random walk would have the vertical coordinate x after time t, 
if there were no motion of the wall. The effect of the moving wall on the winding angle 
distribution is thus a systematic shift in the slopes of the exponential tails. For small values 
of chirality the slopes on the two sides are changed to tt/2 ± a. Due to this explicit velocity 
dependence, these asymmetric distributions are clearly non-universal. At large chiralities, 
the slopes vanish as a//3 2 resulting in quite wide distributions. Apparently strong chirality 
of a defect increases the probability of entanglements. Fig. ^| shows our simulation results 
for the winding angle distribution for a walk on the above mentioned directed triangular 
lattice. The asymmetry due to the shift is clearly visible, and the winding angle distribution 
is wider than for a stationary wall. This case thus exemplifies the strong chirality limit 
discussed in the previous paragraph. We also simulated a square lattice with directed bonds 
as indicated in Fig. 0(b). The corresponding winding angle distribution is shown in Fig. 10. 
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The distribution is again asymmetric, but not as wide as in the previous case, and more 
similar to that expected in the weak chirality regime. 

So far, we have assumed that the winding angle changes by the same amount each time 
the walker returns to the winding center. It is more realistic to assume that the change 
in winding angle has a certain probability distribution. A possible example is provided by 
polymers with randomly changing chirality ||19|| . In the limit t — > oo, the total change in 



winding angle due to the chirality is on an average 

/oo 
d(A9)p(A6)A6, 
-oo 

where p(A6) is the probability distribution for AO. The variance for the total change in the 
scaling variable x is 

and vanishes in the limit t — > oo. The effect of random chirality of the winding angle 
distribution is identical to that of uniform chirality, and is zero when segments of positive 
and negative chirality occur equally often. 

Finally, we want to emphasize that the results of this section are based on the assumption 
that there is no interaction between the polymer and the winding center besides the chirality. 
When there is an additional repulsive interaction, we have to choose absorbing boundary 
conditions, in which case the random walk never hits the moving wall, and the motion of 
the wall (the chirality) has no effect at all on the winding angle distribution. On the other 
hand, when the polymer is bound to the winding center due to an attractive interaction, 
the number of returns to the winding center, and consequently the systematic shift in the 
winding angle due to chirality are proportional to the length t (see section |TTT|) . In the 
presence of both an attractive interaction and a hard-core repulsion (probably the most 
realistic case pOp, the polymer performs a phase transition from a bound to a free state 
depending on the temperature, and we expect the results of this subsection to apply near 
the transition temperature. 
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III. THE WINDING ANGLE DISTRIBUTION OF CONFINED POLYMERS 



Up to now, we considered only cases where the polymer could wander infinitely far 
away from the winding center. However, there are many physical situations where polymers 
indeed are confined to some region around the winding center. When the diameter of the 
container is small compared to the length of a DP, or when the polymer density is so large 
that for each of them just a small cylinder is available, the winding angle distribution will be 
fundamentally different from the previous section. Random walk segements of time At that 
are small compared to the total length of the walk, but large compared to the time it takes 
to make a winding, have identical winding angle distributions. Applying the central limit 
theorem, we can therefore predict that the total winding angle distribution of a confined 
polymer has the form 

(20) 

where 1/a is the variance in the winding angle per unit time. 

When the winding center is chirally asymmetric, there is an additional shift in the mean 
winding angle. In the following three subsections, we will determine the winding angle 
distribution for polymers confined between two cylinders, polymers bound to an attractive 
winding center, and bound polymers winding around chiral centers. 



Pcon(^) oc exp —a9 2 /2t 



A. Polymer confined between two cylinders 

We start with the simple model of a polymer confined between two concentric cylinders 



(see Fig. (TTJ) . The inner cylinder is the winding center, the outer one is the wall of the con- 
tainer, or represents the repulsion of the neighboring polymers. This situation is equivalent 
to a random walk confined between two concentric rings of radii -Ri and i?2- After a long 
time, the probability p(r) to find the walker at a given radius r is independent of time, and 
of the angle. For reflecting boundary conditions at the outer and the inner ring, the walker 



19 



is then with equal probablity at any site between the two rings, leading to 

P(r) = r\n(R 2 /R 1 ) ' 

The variance of the winding angle per unit time is 1/r 2 when the walker is at radius r. Since 
the variances of different segments of the random walk are independent, we can add them 
up, leading in the limit of large t to 

1 f R2 1 f R2 dr t ( 1 1 

-p(r)dr 



a J Rl r 2 ^ v ' J Rl r 3 \n(R 2 /R 1 ) 2\n(R 2 /R 1 ) \R 2 R\) ' 

It is more physically relevant to use absorbing boundary conditions with p(r) = for r = 
Ri,R 2 . We solve the diffusion equation 

dP(r,<j>,t) _ 1 (d 2 P ldP_ }_dP\ 
dt 2 \ dr 2 r dr r 2 <fi 2 J 

with the ansatz 

oo 

P(r, 0, t) = P( r ) cos(n0) exp [-X n t] . 

n=0 

For long times, the angular dependence vanishes, and the mode with the slowest decay (the 
smallest eigenvalue A ) dominates. (Note that is not the winding angle, but the azimutal 
angle, which takes only values between and 2ir.) Since we normalize the winding angle 
distribution with respect to the walkers that do not get absorbed, the factor exp[— \ t] drops 
out, and the probability to find the walker after time t at radius r is given by the solution 
of 

d 2 p(r) 1 dp(r) _ . . 

+ — IT 1 + 2A °P( r ) = °> 

with the boundary conditions given above, and with the normalization condition 

Jri P( r )dr = 1- The general solution of this (Bessel) differential equation can be written in 

form of an integral 

p(r) — J C\ cos ^y / 2Aor sin C^j + C 2 cos ^ \J 2\qt cos C^j In (^J 2X r sin 2 (^j d(. 

The values of Ao, Ci, and C 2 are obtained by matching two consecutive zeros of this function 
to r = Ri and r = R 2 , and by normalizing properly. In general, the solution cannot be 
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written down in a closed form and has to be found numerically. In the case R\ -C -R2, the 
values of Ao and C\jCi can be found analytically, since C\ ^> C2 in this limit, and the first 
two zeros of p(r) are given by the conditions \ia(y/2X Ri) = C\/Ci and a/2Ao-Ri — 2.4 (the 
first zero of the Bessel function J ). In the limit R 2 — Ri <C 1, we findp(r) = C sin (V ) , 
where C is the normalization constant. 

Exercise: Set R\ = 1 and determine numerically the variance 

1 /--R2 p( r )dr 
a J/?! r 2 

as a function of R2 ■ 



B. Polymer bound to an attractive winding center 



A polymer can also be confined by an attractive winding center (see Fig. [121): A DP 
subject to an attractive potential of radius b and binding energy U per unit length, is 
bound to that winding center. For temperatures above the crossover value of T* oc boy/Uo, 
the polymer is only weakly bound and wanders horizontally over a large localization length, 
l±(T) — bo exp[(T/T*) 2 ] ||. The mean vertical distance l z between consecutive intersections 
of the polymer with the defect is consequently proportional to l\. Over this distance, the 
polymer can be approximated by a directed walk which returns to its starting point (the 
winding center) after a time l z . 

Using the result in Eq. ([14]) , we can derive the winding angle distribution p°a{x) for such 
confined random walks. Each walk that returns to its starting point after time t (in the 
z-plane) is composed of two walks of length t/2 going from the starting point to z(t/2). As 
we have seen in subsection P A| , almost all walks of length t/2, when mapped on the plane 
2£/ln(t/2), have their endpoint on a vertical wall at distance one from the origin. The 
winding angle distribution for these walks is given by Eq. (0), with t replaced by t/2. The 
probability that a walk that returns to its starting point has a winding angle 6 is therefore 
obtained by adding the probabilities of all combinations of two walks of length t/2 whose 
winding angles add up to 9, i.e. 
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/oo 
dy 
-oo 



IT 7T 



7T /TO , ,ra, , 

— coth(— )-l , (21) 



4 cosh (717/ ) 4 cosh (tt(x — y)) 

~~ 2sinh 2 (7rx/2) V 2 ™"" v 2 

where x = 20/ ln(t/2i2 2 ). 

For large x ~ x, the above expression decays as x exp [— ttx] . A polymer of length t is 
roughly broken up into t/l z segments between contacts with the attractive columnar defect. 
We can assume that the winding angle of each segment is independently taken from the 
probability distribution in Eq. (|2l| ) with t ~ l z . Adding the winding angle distributions of 
all segments leads to a Gaussian distribution centered around 8 = 0, and with a variance 
proportional to Lln(l z )/l z . 



C. Confined polymers winding around chiral centers 

When the winding center has a chiral asymmetry, the mean winding angle is increased 
by some finite amount AO per unit time. The winding angle distribution is consequently 
modified to 

p c (6) = exp [-5(0 - tA6) 2 /t\ , 

with a mean winding angle proportional to the length of the polymer (see also [pOj] ). 1/5 is 
larger than 1/a, since the variances in the number of returns to the winding center, and in 
AO, both contribute to the variance of the winding angle. For weak chirality, 1/a is close 
to 1/a, and the main effect of the chirality is just a shift of winding angle distribution. For 
strong chirality, we expect 1/a ^> 1/a, and the winding angle distribution becomes very 
broad (similar to the situation discussed in section [11 D| ). 

When the polymer is confined not by a container but by neighboring polymers, it will 
not just wind around one of these neighbors. Kamien and Nelson |22| have shown that when 
chirality is strong, screw dislocations proliferate throughout the polymer crystal. 
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IV. WINDING ANGLES IN RANDOM MEDIA AND FOR SELF- AVOIDING 

POLYMERS 



So far, we only considered winding topologies that can be mapped on ideal random 
walks. However, when the medium in which the polymer is embedded is non homogeneuos, 
the energies of different polymer configurations are different. Examples are polymers in gels 
and porous media [^3] , or magnetic flux lines in high-T c superconductors, pinned by oxygen 
impurities [2T|. We consider the case of quenched randomness, where one end of the polymer 
is fixed 0. 

The behavior of a DP in the presence of short-range correlated randomness is modeled 



by a directed path on a lattice with random bond energies |25| . In 3 or less dimensions, the 
polymer is always pinned at sufficiently long length scales. An important consequence of 
the pinning is that the path wanders away from the origin much more than a random walk, 
its transverse fluctuations scaling as t u , where v « 0.62 in three dimensions, and v = 2/3 in 



two dimensions p6| , |27| . The probability of such paths returning to the winding center are 
thus greatly reduced, and the winding probability distribution is expected to change. 

We examined numerically the windings of a directed path along the diagonal of a cubic 
lattice (see Fig. |5]). To each bond of this lattice was assigned an energy randomly chosen 
between and 1. Since the statistical properties of the pinned path are the same at finite 
and zero temperatures, we determined the winding angle of the path of minimal energy by 



a transfer matrix method. For each realization of randomness, this method |25[] finds the 
minimum energy of all paths terminating at different points, and with different winding 
numbers. This information is then updated from one time step to the next. From each 
realization we thus extract an optimal angle as a function of t. The probability distribution 
is then constructed by examining 2700 different realizations of randomness. To improve the 
statistics, we averaged over positive and negative winding angles. 



The resulting distribution is shown in figure |H| with a scaling variable x = 6/2y\n.t. 
This scaling form is motivated by that of self-avoiding walks, which in two dimensions follow 
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a Gaussian distribution 

Ps 4 = i^) = 7? exp H)- (22) 



The result of the data collapse in Fig. 13 agrees well with the Gaussian distribution 



e \ 11.5 , . _ 2 



Prand \X = T = \ — 6Xp ~1.5x . (23) 

V 2Vlnt/ V 7f v 7 
Directed paths in random media and self-avoiding walks share a number of features which 
make the similarity in their winding angle distributions plausible. Both walks meander 
away with an exponent larger than the random walk value of 1/2. (The exponent of 3/4 for 
self-avoiding walks is larger than v « 0.62 for polymers in 3 dimensions.) As a result, the 
probability of returns to the origin is vanishingly small in the limit t — > oo for both types 
of paths, and the properties of the winding center are expected to be irrelevant. (A simple 
scaling argument suggests that the number of returns to the origin scale as N(t) oc 
The conformal mapping of section [TI| cannot be applied in either case: The density and size 
of impurities in a random medium become coordinate dependent under this mapping, as 
does the excluded volume effect. The winding angle distribution for self- avoiding walks in 
Eq. ( |2"2"D has been calculated using a more sophisticated mapping PJT3|. As a similar exact 



solution is not currently available for polymers in random media, we resort to the scaling 
argument presented next. 

Let us divide the self- avoiding walk, or the directed path, in segments going from t/2 
to t, from i/4 to t/2, etc., down to some cutoff length of the order of the lattice spacing, 
resulting in a total number of segments of the order of hit (see Fig. [14]). The statistical 
self-similarity of the walks suggests that a segment of length t/2 n can be mapped onto a 
segment of length t/2 n+l after rescaling by a factor of 1/2 U . Under this rescaling, the winding 
angle is (statistically speaking) conserved, and consequently all segments have the same 
winding angle distribution. Convoluting the winding angle distributions of all segments, 
and assuming that the correlations between segments do not invalidate the applicability of 
the central limit theorem, leads to a Gaussian distribution with a width proportional to Int. 
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This argument does not work for the random walks considered in section 0, since the finite 
radius of the winding center is a relevant parameter. Different segments of the walk are 
therefore not statistically equivalent, as they see a winding center of different radius after 
rescaling. 

In the pure system, we had to distinguish between repulsive, neutral, and chiral winding 
centers. Since in the presence of point impurities the polymer does not return to the winding 
center as often, these differences are now irrelevant. In fact, it can be shown that even an 
attractive force between the winding center and the polymer cannot bind the polymer to 



the winding center, as long as it does not exceed some finite threshold |[28|| . Perhaps not 
surprisingly, the main conclusion of this section is that the pinning to point randomness 
decreases entanglement. 



V. DISCUSSION AND CONCLUSIONS 

Topological entanglements present strong challenges to our understanding of the dynam- 
ics of polymers and flux lines. In this paper, we examined the windings of a single directed 
polymer around a columnar winding center, or the winding of two DPs around each other. 
By focusing on even this simple physical situation we were able to uncover a variety of in- 
teresting properties: The probability distributions for the winding angles can be classified 
into a number of universality classes characterized by the presence or absence of underlying 
symmetries or boundary conditions. 

For free DPs in a homogeneous medium, we find a number of exponentially decaying 
distributions: If there is no interaction at all between the polymer and the winding center 
(corresponding to reflecting boundaries for random walks) we obtain the distribution in 
Eq. (0). Removing this conservation (absorbing boundaries or repulsive interaction between 
the polymer and the winding center) leads to the distribution in Eq. ( |H| ) whose tails decay 
twice as fast. 

A completely new set of distributions is obtained for chirally asymmetric situations, 
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where the polymer is preferentially twisted in one direction at the winding center. These 
distributions have asymmetric exponential tails, with decay constants that depend on the 
degree of chirality. Strong chirality appears to lead to quite broad distributions. A remaining 
challenge is to find the complete form of this probability distribution by solving the two 
dimensional diffusion equation with moving boundary conditions. 

When the polymer is confined to a finite volume around the winding center, the winding 
angle distribution becomes Gaussian, with a width proportional to the length of the polymer. 
In the presence of chiral asymmetry, the mean winding angle is proportional to the length 
of the polymer, and we have to distinguish again between the limits of weak and strong 
chirality. 

For non-ideal walks, with a vanishing probability to return to the origin, the properties 
of the winding center are expected to be irrelevant. Both self- avoiding walks in d = 2 
dimensions, and polymers pinned by point impurities in d = 3, have wandering exponents 
v larger than 1/2 and fall in this category. We present a scaling argument (supported by 
numerical data) that in this case the probability distribution has a Gaussian form in the 
variable 9/ Vint. Not surprisingly, wandering away from the center reduces entanglement. 
The characteristic width of the Gaussian form is presumably a universal constant that has 
been calculated exactly for self-avoiding walks in d = 2. It would be intersting to see if this 
constant (only estimated numerically for the impurity pinned polymers in d — 3) can be 
related to other universal properties of the walk. Changing the correlations of impurities 
(and hence the exponent v) may provide a way of exploring such dependence. 

There are certainly other universality classes not explored in this paper. For example, we 
did not consider the case of a long-range interaction between the polymer and the winding 
center. Also, the mapping of a DP in a nematic solvent on a random walk is correct only 
to first approximation j|,[5[]. Due to long-range correlations within the nematic solvent, the 
number of returns of the polymer to the winding center does not increase with In t, but with 
ln(lnt) for very large t. 

The results of this paper also provide conjectures for the winding of non directed polymers 
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around a rod. If the self-interaction of the polymer can be neglected for some reason, the 
results of section [TT] can be applied, the parameter r now being the internal coordinate of the 
polymer. To first approximation, this might be correct for a polymer close to the #-point, 
but three-point interactions which ultimately swell the ^-polymer will eventually invalidate 
the result §. 

When the polymer swells to give v > 1/2, as is the case for a self-avoiding random 
walk in three dimensions, the winding center is no longer important. The projection into 
the plane perpendicular to a rod is a walk that wanders away from the winding center 
faster than an ideal random walk. We can therefore apply the results of section |Iy| and 
conclude that the winding angle distribution is a Gaussian, with a variance proportional 
to Int. When the polymer is in the collapsed state, its winding angle distribution is again 
different, and has still to be found. Since collapsed polymers are relatively compact, they 
can be approximated by Hamiltonian walks that visit each site within a volume of the size 
of the polymer exactly once |29],|30|. We thus conclude this review with the open problem 



of determining the probability distribution for windings in the collapsed state by examining 
the behavior of Hamiltonian walks. 
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FIG. 1. a) Two directed polymers winding around each other, b) A directed polymer winding 
around a rod. 

FIG. 2. Conformal mapping of the random walk with absorbing boundary conditions, and sub- 
sequent rescaling. 
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FIG. 3. Winding angle distribution for random walks on a square lattice with reflecting boundary 
conditions for t = 38 (dot-dashed), 152 (long dashed), 608 (dashed), 2432 (dotted), and 9728 (solid). 
The horizontal axis is x = 2#/ln(2i). The thick solid line is the analytical result of Eq. (|j~7|) . 

FIG. 4. Winding angle distribution for random walks on a square lattice with absorbing boundary 
conditions. The symbols and the variable x are the same as in the previous figure. The thick solid 
line is the analytical result of Eq. (|i"4| ) . 

FIG. 5. a) Construction of a polymer directed along the (l,l,l)-diagonal of a cubic lattice, b) 
projection into the plane perpendicular to the (l,l,l)-diagonal. The three sublattices are indicated 
by different shades of gray. 

FIG. 6. Winding angle distribution around the preferred direction for a flux line (directed path) 
in 3 dimensions for t = 243 (dot-dashed), 729 (long dashed), 2187 (dashed), 6561 (dotted), and 
19684 (solid). The horizontal axis is x = 2(9/ln(2t). The thick solid line is the analytical result of 
Eq. ©. 

FIG. 7. Triangular and square lattices with directed bonds. The winding centers are indicated 
by a circle o. 

FIG. 8. Brownian motion around a rotating winding center. 

FIG. 9. Winding angle distribution around the preferred direction for a random walk on a di- 
rected triangular lattice for t = 243 (dot-dashed), 729 (long dashed), 2187 (dashed), 6561 (dotted), 
and 19684 (solid). The scaling variable is x = 2#/ln(2i). The thick solid line is the distribution 
given in Eq. (|l7|). 

FIG. 10. Winding angle distribution for a random walk on a directed square lattice for t = 38 
(dot-dashed), 152 (long dashed), 608 (dashed), 2432 (dotted), and 9728 (solid). The thick solid 
line is the distribution given in Eq. (|i~7|). 
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FIG. 11. A directed polymer confined between two cylinders. 

FIG. 12. A directed polymer bound to an attractive rod. 

FIG. 13. Winding angle distribution for a directed path in a random 3-dimensional system for 
t = 120 (dotted), 240 (dashed), 480 (long dashed), 960 (dot-dashed), and 1920 (solid). The thick 
solid line is the Gaussian distribution in Eq. (p3|) . 

FIG. 14. Division of the self-avoiding walk into self-similar segments. 
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